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Elastic  Constants  of  Quartz 


I.  INTRODUCTION 

The  critical  performance  characteristics  of  a  crystal  resonator  device,  fre¬ 
quency.  electromechanical  coupling,  and  temperature  sensitivity,  are  determined 
by  the  elastic,  piezoelectric,  and  dielectric  constants  of  the  material  and  their 
temperature  coefficients.  In  this  investigation  we  consider  the  influence  of  the 
material  coefficients  on  the  thickness  vibrations  of  crystallographically  doubly 
rotated  quartz.  The  simplest  mathematical  formalism  that  describes  a  vibrating 
infinite  plate  with  appropriate  boundary  conditions  relates  the  measured  fre¬ 
quencies  to  the  material  constants.  The  set  of  four  coupled  differential  equations 
yields  a  set  of  linear  equations  for  the  displacement  and  results  in  three  eigen¬ 
values  for  the  phase  velocities.  The  general  formulation  of  this  approach  is  out¬ 
lined,  for  example,  in  Tiersten.  * 

For  an  arbitrary  doubly  rotated  quartz  cut,  the  vibration  frequencies  are 
determined  by  ten  independent  fundamental  material  constants:  six  elastic,  two 
piezoelectric,  and  two  dielectric.  In  turn,  the  vibration  frequencies  can  be  used 
in  conjunction  with  the  mathematical  formalism  to  derive  numerical  values  for  the 

(Received  for  publication  13  May  1982) 

1.  Tiersten,  H.F.  (1969)  Linear  Piezoelectric  Plate  Vibrations,  Plenum  Press, 
pp.  88-93. 


material  properties.  The  derived  material  values  are  mathematical-formalism 
dependent,  and  they  may  differ  from  the  intrinsic  constants. 

For  highly  symmetric  crystal  orientations,  the  cubic  equation  for  the  fre¬ 
quency  eigenvalues  can  be  factored,  and  the  number  of  material  constants 
involved  in  determining  the  vibration  frequencies  is  drastically  reduced.  In  some 
cases  the  resulting  frequency  depends  on  a  single  elastic  constant.  These  direc¬ 
tions.  designated  as  the  uncoupled  mode  orientations,  are  the  ones  most  often 
utilized  for  the  experimental  determination  of  the  elastic  constants. 

The  experimental  determination  of  the  elastic  constants  of  quartz  has  occupied 

2 

many  investigators  over  long  periods  of  time.  Cady  reviews  several  approaches 
and  lists  corresponding  numerical  values.  Reference  2  is  a  corrected  reprint  of 
the  original  treatise  published  on  the  subject  of  piezoelectricity  in  1946,  and  the 

listed  values  are  the  ones  accepted  in  1946.  Currently,  the  most  widely  accepted 

3 

elastic  constant  values  are  the  ones  determined  by  McSkimin,  and  independently  by 
Bechmann,  Ballato.  and  Lukaszek.  ^  [This  reference  is  hereafter  abbreviated  as 
BBL  (1962)]  .  McSkimin  has  applied  ultrasonic  excitation  to  bulk  crystals  along 
the  uncoupled  directions,  whereas  BBL  utilized  piezoelectric  excitation  and  doubly 
rotated  crystal  disks  in  a  resonator  configuration.  Both  investigators  used  the 
same  mathematical  description  to  relate  experimental  wave  velocities  to  the 
elastic  constants.  The  elastic  constants  determined  by  these  investigators  agree 
with  each  other,  and  also  generally  agree  with  previous  results. 

We  have  applied  a  least  squares  analysis  to  the  BBL  (1962)  data  set  and  we 
find  major  discrepancies  between  the  listed  values  and  the  least  squares  fit  derived 
constants.  This  raises  serious  questions  with  respect  to  the  validity  of  the  mate¬ 
rial  values  and  their  applicability  to  crystallographically  doubly  rotated,  piezo- 
eleetrically  excited  resonators. 


2.  THICKNESS  VIBRATIONS  OF  PIEZOELECTRIC  PLATES 

Tiersten1  has  derived  the  differential  equations  for  the  thickness  vibrations  of 
an  arbitrarily  anisotropic,  piezoelectric,  infinite  plate.  The  wave  propagation  is 
described  by  a  set  of  four  coupled  linear  equations  in  four  variables  u.  and  ® 


2.  Cady.  W.G.  (1964)  Piezoelectricity.  Dover  Publications,  Vol.  1,  pp.  134-157. 

3.  McSkimin,  H.  J.  (1962)  Measurement  of  the  25°C  zero-field  electric  moduli 

of  quartz  by  high  frequency  plane-wave  propagation,  J.  Acoust.  Soc.  Am. 
34:1271-1274.  '  . . 

4.  Bechmann,  R.  ,  Ballato,  A.D. ,  and  Lukaszek,  T.J.  (1962)  Higher-order 

temperature  coefficients  of  the  elastic  stiffnesses  and  compliances  of 
alpha-quartz,  Proc.  of  the  IRE  50:1812-1822. 
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where  cj^k(/,  e'VVj.  and  £},„  are  ihe  elastic,  piezoelectric,  and  dielectric  con¬ 
stants,  p  denotes  the  density,  vl  the  components  of  the  displacement  vector,  *  the 
electromagnetic  scalar  potential,  xv  the  spatial  coordinates,  and  t  the  time.  The 
primed  notation  on  the  material  constants  indicate  that  these  quantities  are  referred 
to  the  rotated  wave  propagation  coordinate  system.  A  standard  tensor  convention 
is  used  in  subscripts,  namely,  repeated  Latin  indices  imply  summation  over  1,  2, 
and  3,  but  no  summation  is  implied  over  repeated  Greek  indices.  In  succeeding 
equations  we  also  adopt  the  notations  of  indicating  differentiation  with  respect  to 
space  coordinates  by  commas  followed  by  an  index  and  time  derivatives  by  dots 
over  the  variables. 

The  boundary  conditions  applicable  to  this  problem,  a  forcing  alternating 
voltage  and  vanishing  mechanical  stresses  at  the  plate  surfaces,  x^  =  ±h,  are 


(3) 

(4) 

where  T.  is  the  stress  tensor,  and  2h,  the  plate  thickness. 

The  piezoelectrically  "stiffened”  elastic  constant  is  defined 


®  =  ±  ®  cos  ut 
o 


T  =  c'  _ !i  +  e'  d*  =  0 

Ki  9x^  m/j  ax^ 


icjk  r 


-  C  *  -f  P*  p 

jcjk  v  vvj  i/vk 


/*{ 


(5) 


In  terms  of  „, 
equations 


Eqs.  (1)  and  (2)  combine  to  form  a  set  of  three  linear  coupled 


v 


pUi 


(6) 


Tiersten  shows  that  a  solution  to  the  differential  equations  based  on  a  linear  com¬ 
bination  of  symmetric  and  antisymmetric  plain  waves  remains  uncoupled,  and, 
moreover,  that  the  symmetric  waves  cannot  be  excited  by  the  application  of  an 
alternating  voltage  to  the  surface  electrodes.  Consequently,  it  is  mathematically 
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rigorous  to  assume  that  the  solution  is  a  set  of  antisymmetric  plain  waves  in  the 
form 


u^Xy,  t)  =  UjiXy)  cos  ut  =  sin  (w/v)  x^  cos  ut  (7) 

Kx^,  t)  =  Q(*v)  cos  <0*.  (8) 

Substituting  these  expressions  into  Eq.  (6),  yields  the  set  of  three  equations 

«Vjkv-6jkr)0k  =  O  3.  ^  '  1,2,3  (9) 

where  the  eigenvalue  r  is  related  to  the  phase  velocity  v  by 

r  =  pv2  (10) 


The  three  eigenvalues  r  are  obtained  by  expanding  the  3X3  determinant 


'v,  -  6jkr! =  0 


an 


and  solving  the  resulting  cubic  equation  for  r  m,  m  =  1,  2,  and  3,  For  a  given 
frequency  u  there  are  three  independent  transverse  velocities. 

For  a  given  eigenvalue  r™  and  eigenvector  (1™,  the  electromagnetic  potential 
is  determined  from  the  boundary  conditions  of  Eq.  (3).  This  results  in 


where  the  integration  constant  Lj  is 

4  =  *0  -  fT  sin  (uh/vm)  (13> 

The  short-circuit  resonance  frequencies  are  found  from  the  six  stress  boundary 
conditions  implied  in  Eq.  (4).  They  are  the  solutions  of  the  transcendental  equation 


m  - 1 


tan  (wh/vm) 
<toh/vm) 


(14) 


whcr.-  for  a  given  mode  m  the  electromechanical  coupling  factor  km  is  defined  as 
(km>2  =  <e^  ^>2/(eUrm)  m  =  1.2.3  (15) 
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The  sum  e'  .  dm  is  the  effective  piezoelectric  constant,  and  it  determines  the 

t't/j  ,1 

excitation  magnitude  of  the  mode.  In  general,  the  resonance  frequencies  cannot 
be  attr  ibuted  to  a  particular  mode,  but  are  coupled  functions  of  all  three  modes. 
Ballato3  reviewed  the  regions  in  which  different  approximations  applicable  to  the 
resonance  equation  are  valid. 


3.  COORDINATE.  SYSTEMS  AND  TRANSFORMATION  MATRICES 


Figure  1  depicts  the  plate  geometry  and  crystallographic  coordinate  system 
used  in  this  investigation.  We  assume  Xj  =  x,  Xg  =  v,  and  x^  =  z,  with  Xj  and 
Xy  in  the  plane  of  the  plate.  The  phase  velocity  propagates  along  Xj,  normal  to 
the  plate.  The  initial  position  of  the  plate  is  a  Y-cut,  w  ith  its  thickness  direction 
normal  to  the  y-axis.  The  position  of  the  rotated  plate  is  described  by  angles 
(<t>,  0).  where  <f>  and  0  are  rotations  around  the  z-  and  x-axis,  respectively. 


Figure  1.  Coordinate  System  for  a  Doubly  Rotated  Plate 


The  wave  propagation  and  crystallographic  coordinate  system  3  are  related  by 
a  transformation  matrix  V....  The  matrix  elements  are  defined  by  the  dot 
products  of  the  orthogonal  unit  vectors  of  the  two  systems.  The  direction  cosines 
for  a  ^-rotation  around  the  z-axis  are  given  by 


aii  = 


cos  <j>  sin  <t>  0 

-sin  <f>  cos  d  0 

0  0  1 


(16) 


5.  Ballato,  A.  (1977)  Doubly  rotated  thickness  mode  plate  vibrators.  Physical 

Acoustics.  XIII,  W.  P.  Mason  and  R.  N.  Thurston,  Eds. ,  Academic  Press, 
New  York,  pp.  115-181. 
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and  for  a  0-rotation  around  the  x-axis  by 


1  0  0 

0  cos  0  sin  0 

0  -sin  8  cos  6 


(17) 


For  a  doubly  rotated  crystal  plate,  a  ^ -rotation  followed  by  a  0-rotation,  the 
transformation  matrix  V . ^  =  bj^a^j  connecting  the  wave  propagation  and  crystal¬ 
lographic  coordinate  systems  becomes 


V;. 


cos  <j> 

-sin  $  cos  0 
sin  <j>  sin  0 


sin  ^  0 

cos  ij>  cos  0  sin  0 
-cos  sin  0  cos  0 


(18) 


The  material  properties  appearing  in  the  differential  equations  are  described 
by  2,  3,  or  4  subscripts.  We  introduce  the  compressed  subscript  notation,  replac¬ 
ing  ii  or  kl  by  p  or  q,  according  to  the  prescription  given  in  Table  1.  This  allows 
the  material  constants  to  be  described  by  two  subscripts  and  expressed  as 
mat  rices. 


Table  1.  Relationship  Between 
Subscripts  i.  j,  k,  1,  and  p,  q 


The  primed  material  constants  referred  to  in  the  wave  propagation  system  are 
related  to  quantities  in  the  crystallographic  reference  system  by  second,  third, 
or  fourth  rank  tensor  transformations.  More  specifically. 


V. 

im 


V.  e 
.in'  mn 


(19) 
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(20) 


? 

! 


e!  =  e!.,  =  V.  V.  V.  e 
ip  ljk  im  jn  ko  mno 

c'=c'  =V  VVVc 
pq  ijkl  im  jn  ko  lp  mnop 

where  V-  are  the  rotation  matrix  elements  of  Eq.  (18),  and  the  repeated  sub¬ 
script  indices  imply  summations  over  1,  2,  and  3.  Thus  each  fU,  ej  ,  and  cj^ 
element  is  the  sum  of  9,  27,  and  81  terms,  respectively. 

In  the  pq  notation,  e  is  described  by  a  3  X  3  matrix  with  9  elements,  e  by  a 
3X6  matrix  with  18  elements,  and  c  by  a  6  X  6  matrix  with  36  elements.  How¬ 
ever,  for  any  crystal  system  the  number  of  independent  matrix  elements  is  re¬ 
duced  by  symmetry.  For  example,  for  the  most  general  crystal  system,  a 
triclinic  system  without  a  center  of  symmetry,  the  applicable  symmetries  are: 


£ij  £ji 


Pijk  =  Pikj 


cijkl  °ijlk  cjikl  °klij  °f  1  pq  ‘  qp 


(22) 

(23) 

(24) 


The  number  of  independent  constants  arc  thus  reduced  to  45  (6  dielectric,  18 
piezoelectric,  and  21  elastic).  The  number  of  independent  material  constants  are 
further  reduced  with  increasing  crystal  symmetry. 

The  material  of  primary  interest  to  this  investigation  iso -quartz.  This 
trigonal  crystal  system,  (32,  class  D^),  has  a  three-fold  symmetry  with  respect 
to  the  z-axis  and  a  two-fold  symmetry  with  respect  to  the  x-axis.  The  material 
properties  are  invariant  under  a  120°  rotation  around  the  z-axis  or  a  180°  rotation 
around  the  x-axis.  Using  these  transformation  properties,  it  can  be  shown  that  the 
number  of  independent  material  constants  reduce  to  10,  2  dielectric,  2  piezoelec¬ 
tric,  and  6  elastic.  The  material  matrix  arrays  applicable  too-quartz  are  listed 
in  Table  2.  The  structure  of  these  arrays  also  holds  true  for  other  material 
properties  transforming  as  second,  third,  or  fourth  order  tensors.  For  example, 
the  array  for  the  linear  expansion  coefficient  o,,  has  the  same  matrix  structure  as 


Table  2.  Material  Arrays  for  o  -Quartz 


One  can  derive  explicit  algebraic  expressions  for  the  "rotated”  material  con¬ 
stants  of  Eqs.  (19),  (20),  and  (21),  by  applying  the  rotation  matrix  V.j  to  the  arrays 
listed  in  Table  2.  The  results,  for  the  material  components  of  interest  to  this 
investigation,  are  listed  in  Table  3.  For  quartz,  owing  to  crystal  symmetry,  all 
rotations  can  be  mapped  into 

0  <  *  s  60°  and  0  <  0  <  90° 

In  this  nomenclature,  the  standard  AT-cut  crystal  is  designated  as  (0,  ~35)°  and  the 
BT-cut  as  (fi0,~49)°.  Another  commonly  used  mapping  region  is 


14 


and 


0  ■£  <t>  s  30 


-90  <  9  <  90u 


The  two  descriptions  are  related  by  the  transformations 
-  60°  -  4  and  9  -»  -0 

In  this  nomenclature  the  BT-cut  becomes  (0,  -49)°,  and,  for  example,  the  com¬ 
bination  (40° 54',  16°34')  transforms  to  (19°0G',  -16°34'). 


Table  3.  Relationship  Between  "Rotated"  and  Crystal -Axes -Referred  Material 
Constants  for  Quartz 


2  2 
e  j  j  cos  9  +  sin  9 


-  e^j  sin  3  <j>  cos'  0 

2 

-e.j  sin  3  <t>  sin  9  cos  9 
=  -(e1  j  cos  3  4  cos  9  +  e14  sin  0)  cos  6 

-  Cjj  cos4  9  +  c.j.j  sin4  9  +  j 3 / 2  +  e^)  sin* 2  29 
-  2  c  l4  cos  3  4  sin  2  9  cos2  9 

(l/alKcj^  +  c33>  sin2  9  +  (c13  +  2  c44>  cos  2  9  -  c^J  sin  2  9 

'  r 14  ros  ^ ^  cos  ^  cos  9  9 

2 

=  -3  Cj4  sin  3 d>  sin9  cos  9 

=  c44  +  <cll^4  ■  c  13^2  +  c33/4  -  c44>  sl°2  2  9  +  <l/2>  C14  cos  3^  sin  49 
=  cH  (2-3  cos2  9)  sin  3  4  cos  9 
2  2 

=  cfifi  cos  9  +  c44  sin  9  +  cJ4  cos  3  4  sin  2  9 
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4.  THICKNESS  VIBRATIONS  OF  QUARTZ 


Consider  thickness  vibrations  normal  to  the  Xg-axis.  The  differential  equa¬ 
tions  of  motions,  Eq.  (G),  expand  to 

rG6  ul,  22  +  r26  u2,  22  +  F46  u3,  22  =  pul 

r26  ul.  22  +  r22  u2,  22  +  F24  u3,  22  =  pu2  (2: 

r46  ul.  22  +  r24  u2,  22  +  r44  u3,  22  =  pu3 

The  appropriate  "stiffened"  elastic  constants,  referred  to  the  crystal  axes,  are 
listed  in  Table  4.  The  three  eigenvalues  Tm  are  calculated  by  expanding  the 
determinant  derived  from  Eq.  (11). 

Table  4.  "Stiffened"  Elastic  Constants  of  Quartz 

r22  =  fcn  +  (e^/e)  sin2  3  ^1  cos4  0  +  Cj3  sin4  0  +  (cJ3/2  +  c^4>  sin2  2  0 
-  2  Cj4  sin  2  0  cos2  0  cos  3^ 

r44  c44  cos2  2  0  +  (1/4MCJ  j  +  (e^/e)  sin2  -2c[3  +  c^]  sin2  20 

+  (1/2)  c14  sin  40  cos  3^ 

+  (e2j/e)  cos2  3  <»j  cos2  0  +  (c^  +  e24/e)  sin2  0 
+  (c^4  +  Cjj  c14/e)  sin  20  cos  3$ 
r24  ({-(cjj  +  fejj/e)  sin2  3^J  cos2  0  (c13  +  2  c44)  cos  2  0 

+  c33  sin2  0}  sin0  -  c  14  cos  3  4  cos  3  0)  cos  0 
r2G  =  -  1(3  c  j4  +  ej  1  e14/e)  sin  0  +  (e^/e)  cos  3  tj>  cos  0]  sin  3  $  cos2  0 
r4fi  =  ((3  sin2  0  -  1)  c14  +  (l/2)(e21/e)  cos  3  ^  sin  2  0  +  (e^  eI4/e)  sin2  01 
sin  3  tj>  cos  0 

where 

«  =  *11  +  e33  tan*® 
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’I 


pq 


into 


r66  '  r 

r26 

r46 

r26 

r22  ‘  r 

r24 

=  0  , 

(26) 

r46 

r24 

r44’r 

!  +  A^  r  + 

o 

II 

o 

(27) 

where  - A ^  is  the  trace,  A^  the  sum  of  the  minors  of  the  major  diagonals,  and 
-Aq  the  value  of  the  determinant.  The  coefficients  are  given  by 


A2  =  '(r22  +  r44  + 

A1  =  (r22r44  "  r24)  +  (r22r66  *  r>26)  +  <r44r66  *  r46) 

Ao  =  r22r46  +  r44r26  +  r66r24  *  r22r44r66  *  2r24r26r46 


(28) 

(29) 

(30) 


Equation  (27)  is  solved  for  the  three  eigenvalues  rm  by  standard  algebraic  proce¬ 
dures  applicable  to  a  cubic  equation  with  three  real  roots. 

The  three  numerically  ordered  phase  velocities  are  usually  designated  as 


v  >  v.  >  v 
a  b  c 


and  are  calculated  from  Eq.  (10).  Figures  2a  through  2g  show  the  angular  varia¬ 
tions  of  wave  velocities  for  $  =  0  to  60°  in  intervals  of  10°.  For  these  calculations 
we  have  used  the  c^,  e.^,  and  values  of  BBL  (1962).  The  types  of  plots  shown 
in  Figures  2a  through  2g  were  initially  published  by  Bechmann  in  1935,  republished 
in  Kef.  6,  and  reprinted  in  Ref.  5.  We  include  them  here  for  reference  and  illus¬ 
tration.  The  numerically  ordered  mode  designations  v  ,  v^,  and  v  are  fairly 
arbitrary,  as  for  a  given  <t>  the  velocity  curves  may  cross  as  a  function  of  0.  This 
can  be  seen,  for  example,  for  d  =  60°,  where  the  b-  and  c-modes  interchange  at 
0  —  25°.  The  right  hand  side  scales  of  Figures  2a  through  2g  are  the  frequency- 
thickness  constants,  defined  as  vm/2,  and  they  can  be  used  to  calculate  an  approx¬ 
imate  plate  thickness  for  a  specified  fundamental  or  overtone  frequency. 

6.  Bechmann,  R.,  Ballato,  A.D.,  and  Lukaszek,  T.J,  (1961)  Frequency- 

temperature  behavior  of  thickness  modes  of  doubly-rotated  quartz  plates, 

15th  Annual  Symposium  on  Frequency  Control,  AD  265455. 
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ANGlf  Of  ROTATION.  8IDTG> 


igure  2a.  Propagation  Velocity,  Frequency  Constant,  and  Cou 
ing  Factor  for  Doubly  Rotated  Quartz,  0 


A\Gtl  Of  RO1ATI0N,  »IDIG' 


Figure  2b.  Propagation  Velocity,  Frequency  Constant,  and  Cou 
pling  Factor  for  Doubly  Rotated  Quartz.  4  10J 


ANCLE  OF  ROTATION,  •  (DEG) 


Figure  2e.  Propagation  Velocity,  Frequency  Constant,  and  Cou 
pling  Factor  for  Doubly  Rotated  Quartz,  $  -  40° 


Figure  2f.  Propagation  Velocity,  Frequency  Constant,  and  Cou¬ 
pling  Factor  for  Doubly  Rotated  Quartz,  ^  =  50° 


ANGIE  OF  ROTATION.  MDEGI 


Figure  2g.  Propagation  Velocity,  Frequency  Constant,  and  Cou¬ 
pling  Factor  for  Doubly  Rotated  Quartz,  <j>  =  60 


The  calculation  of  the  electromechanical  coupling  factor  [Eq.  (15)]  includes 

both  the  eigenvalues  and  the  eigenvectors.  Other  eigenvalue-  and  eigenvector- 

dependent  field  quantities  include  stress,  strain,  electric  displacement,  and  power 

density.  Explicit  formulas  for  the  field  quantities  are  given  by  Slobodnik, 

7 

Delmonico,  and  Conway.  Standard  mathematical  subroutines  are  available  at 
most  computer  installations  for  calculating  eigenvalues  and  eigenvectors  of  any 
matrix  type  and  size.  For  the  3X3  real,  symmetric  matrix  applicable  to  this 
problem,  one  can  derive  analytic  expressions  for  the  eigenvalues  and  eigenvectors. 
This  simplifies  programming  requirements,  and  computations  can  be  performed 
even  on  hand-held  or  desk-top  calculators. 

The  eigenvectors,  satisfying  the  orthonormality 


|3"  P™  =  6 
'11  nm 


are  derived  by  expanding  the  determinant  |r  |  cofactors  and  constructing  the 
adjoint, 


Slobodnik,  A.J.,  Delmonico,  R.T.,  and  Conway,  E.  D.  (1980)  Microwave 
Acoustic  Handbook,  Vol.  3.  Bulk  Wave  Velocities,  RADC-TR-80-188, 
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The  equations  derived  in  this  section  are  analytic  and  self-consistent.  Given 
the  material  constants,  one  calculates  for  the  particular  angular  orientation  (tj>,  6) 
of  interest,  the  "stiffened"  elastic  constants  from  Table  4,  solves  the  cubic  equa¬ 
tion,  Eq.  (27),  for  the  eigenvalues  rm,  and  finds  the  coupling  factors  km  from 
Eqs.  (33)-(35).  This  then  provides  a  convenient  formalism  to  calculate  the  thick¬ 
ness  vibration  propagation  velocities  and  electromechanical  coupling  factors  for 
any  doubly  rotated  quartz  plate. 

In  this  section  we  have  derived  the  velocities  by  considering  the  wave  propa¬ 
gation  along  the  x„-axis.  An  alternate  formulation,  carried  out  by  Slobodnik 
7  i 

et  al,  considers  the  transverse  propagation  along  the  Xj-axis.  The  results  are 
identical,  but  the  pq  subscripts  differ.  The  relationship  between  the  two  formula¬ 
tions  is  given  in  Appendix  A. 


5.  UNCOUPLED  MODES  OF  VIBRATION 


The  eigenvalue  determinant  [Eq.  (26)]  can  also  be  expanded  into 


<r22  *  r)(r44  -  r><r6B  '  n  -  r46<r22  ’  r)  ‘  r26(r44  '  F) 


r24(r66  "  r)  '  2  r24r26r4li  ""  0 


(36) 


When  any  two  of  the  three  off-diagonal  elements  vanish,  the  equation  factors  into 
a  linear  and  quadr  atic  term  in  T.  The  factored  equations  are  listed  in  Table  5. 

If  two  of  the  off-diagonal  elements  vanish,  the  three  vibration  eigenvalues 
corresponding  to  m  -  1,  2,  and  3,  are  given  by 


r(1):  r(2)  =  U/2)<rpp  +  rqq>  ±  <i/2>i<rpp  -  rqq>2  +  4r2q)  1/2 


(37) 


and 


,(3) 


(38) 


where  p,  q,  and  r  ar  e  permutations  of  2,  4,  and  6.  The  (<t>,  6)  combinations  that 
satisfy  Tpq  -  0  can  be  evaluated  from  the  expressions  given  in  Table  4.  Figure  3 
is  a  plot  of  these  solutions  using  the  material  constants  given  in  BBL  (1962).  In 
Figure  3,  the  overlapping  curves  at  4  =  0  and  GO0  and  at  6  =  90°  are  displaced  for 
clarity. 


Table  b.  Factors  of  Wave  Velocity  Eigenvalue  Equation 


90  ■ 

75 

60  ■ 

45 

30  • 

15  ■ 

0  ■ 

0  20  40  60 


4> 


=  0 
=  0 
=  0 


Figure  3.  Zero  Values  of  Off-Diagonal 
Matrix  Elements,  r  =0 

pq 


We  note  that  T24  =  0  is  satisfied  by  (d,  90)°,  that  is,  for  any  <t>  at  9  =  90°.  It 
is  also  satisfied  by  a  continuous  function  between  <t>  -  0  and  60°,  with  end  points 
(0,  40.  8)°  and  (60.  72.  4)°,  and  by  a  third  branch  continuous  between  (30,  0)°  and 
(60,  17.  7)°. 

r20  r  0  is  satisfied  by  any  6  at  4  =  0  and  60°,  and  for  any  <t>  at  0  =  90°,  that  is, 
(0,  0)°,  (60,  9)°.  and  (6.  90)°.  In  addition,  the  bracketed  term  of  r„„,  listed  in 
Table  4,  yields  a  slowly  varying  function,  symmetric  around  (30,  0)°  and  increas¬ 
ing  to  6  0.  8°  at  <t>  =  0  and  60°.  The  deviation  from  (£,  0)°  is  due  to  the  small 

contribution  of  the  piezoelectric  constants  to  the  bracketed  term. 
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r46  =  0  has  solutions  (0,  0)°.  (60,  0)°,  and  (d.  90)°.  The  bracketed  term 
yields  a  continuous  function  between  (0,  35.  5)°  and  (60,  34.  9)°.  The  small  varia¬ 
tions  in  0  are  again  due  to  piezoelectric  effects.  If  we  neglect  piezoelectric 
effects  we  obtain  (d,  35.  26)°,  which  is  exactly  the  AT-cut  angle.  It  is  uncertain 
whether  this  is  a  pure  coincidence,  or  is  a  contributing  factor  in  determining  the 
position  of  the  AT-cut. 

The  overlap,  or  intersection,  of  two  or  three  solutions  of  r  =  0  gives  the 
uncoupled  modes  of  vibration.  Prom  the  diagram  depicted  in  Figure  3,  these 
directions  are  given  by  (0,  0)°,  (60,  6)°,  (<t>,  90)°,  and  (30,  0)°.  One  of  the  problems 
associated  with  the  preferred  directions  is  that  the  piezoelectrically  "stiffened" 
elastic  constant,  defined  by  Eq.  (5),  may  no  longer  be  "stiffened",  has  no  piezo¬ 
electric  coupling,  and  cannot  be  excited  by  piezoelectric  means.  An  inspection  of 
Table  4  shows  that  for  and  the  term  containing  the  piezoelectric 

constant  is  multipled  by  a  sin^  3  d  factor,  which  vanishes  both  at  d  =  0  and  60°. 
Thus  at  (0,  8)°  and  (60,  6)°  two  of  the  modes  are  piezoelectrically  inactive,  and  at 
(d.  90)°  none  of  the  modes  can  be  excited  by  an  alternating  voltage  applied  to  the 
electrodes, 

3 

McSkimin  has  evaluated  the  elastic  constants  from  ultrasonic  excitation. 

He  selected  (0,  0)°,  (0,  90)°,  and  (30,  0)°.  At  these  points,  the  three  modes  of 
vibrations  reduce  to  values  listed  in  Table  6.  The  choice  of  (0,  90)°  is  not  unique. 
Any  (d,  90)°  will  give  a  doubly  degenerate  3  c. .  and  T ^  =  c„„. 

O  O 

Similarly,  (0.0)  and  (60,  0)  yield  identical  values.  From  the  angular  orientation 

2 

chosen  by  McSkimin  one  can  evaluate  e^/e^  add  all  elastic  constants  except  c^g. 
The  determination  of  c^  involves  four  other  elastic  constants,  and,  depending  on 
accumulation  or  cancellation  of  experimental  errors,  may  be  the  least  accurate. 

By  the  same  token,  Cgg  and  c44  should  be  very  accurate.  For  Cj3  McSkimin 
chose  (0,  45)°,  at  which  point 

c13  ^c14  "  c44^  +  “  ^c33  +  c44>'[2r  *  ^  c14  +  c44^1  l  ^  (39) 

but  any  other  angular  orientation  will  result  in  similar  complexity. 

If  one  is  restricted  to  determining  the  material  parameters  from  piezoelectric 
excitation,  for  example,  from  frequencies  derived  from  resonators  or  disks  in 
dielectric  gap  holders,  the  key  equation  determining  the  angular  cuts  is 


r(3)<ro.0>;  rO)(60,  e )  --  <c66  +  e2n/e)  cos2  e  +  (c44  +  <^4/e>  sin2  0 


where  the  plus  sign  applies  to  6  =  0,  and  the  minus  sign  to  6  =  60 
selection  of  6  at  6  n 


The  proper 


0  and  60u  will  yield  c14,  c44,  c^g,  the  piezoelectric,  and  the 
dielectric  constants,  it  is  especially  important  to  have  identical  6-angle  cuts  for 
both  6-0  and  60°.  This  allows  simplifications  in  terms  of  r^(0,  9)  ±r*3^(60,  9). 
As  before,  Cjj  is  determined  from  (30,  0)°,  but  neither  Cj^  nor  c^g  are  accessible 
to  uncoupled  mode  piezoelectric  excitation. 


Table  6.  Eigenvalues  rm  for  Highly  Symmetric  (6,  9)  Orientations 


(6.  9) 

!  rfl>;  r<2) 

;  r(3> 

(0.  0)  and  (60,  0) 

!  (l/2)((cu  +  c44)±[(cu-c44)2  +  4c24 

1 

1/2 ! 

1  ^  2  / 

1  c66  +  ell/£ll 

1 

(6,  90) 

i  C'44 

]  c33 

(30.  0) 

|  (l/2)((c44  +  cRR)  ±((c44-c(.(.)2  +  4c“4 

1 

|l/2} 

1  2 

1  ‘'ii  +  en/en 

1 

6.  EXPERIMENTAL  VALUES  OF  ELASTIC  CONSTANTS  OF  QUARTZ 

Cur  lent  lv,  the  most  widely  accepted  c  _  values  are  the  ones  derived  by 
•i  '  P9 

McSkimin  and  independently  by  BBL  (1962),  and  these  are  in  general  agreement 
with  previous  investigations.  In  Table  7,  under  the  heading  "Observed",  we  list 
the  experimental  data  sets  of  both  investigations,  and  in  Table  8,  we  list  the 
derived  elastic  coefficients.  The  differences  between  the  McSkimin  and  BBL  de¬ 
rived  values  are  small,  and  either  set  can  be  used  equally  well  to  calculate  appro¬ 
priate  quantities  of  interest.  In  Table  8,  we  also  include  the  other  material 
properties  needed  for  calculating  the  transverse  wave  propagation  of  quartz.  Unless 
otherwise  specified,  these  are  taken  from  BBL  (1962). 

In  Table  7,  under  the  heading  "BBL  (1962)”.  we  list  three  sets  of  calculated 
values.  The  first  column  is  taken  from  BBL  (1962),  and  the  second  column, 

"This  Work  1",  is  calculated  in  this  investigation,  using  BBL  (1962)  derived  con¬ 
stants.  The  slight  differences  in  values  are  due  to  the  fact  that  BBL  (1962)  calcula¬ 
tions  are  based  on  a  "semi-stiffened"  elastic  constant, 

r  -  ((■  +  e„  p.  /f,,)1 

pq  pq  2p  l  q'  2  Z 
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rather  than  the  rigorous  definition,  Eq.  (5) 


t 


pq 


c'  + 

pq 


e2pe2q/e22 


Table  7.  Measured  and  Calculated  Velocities  for  Doubly  Rotated  Quartz 
Orientations 


Orientation 
Angles 
<t>  0 

I  Mode 

1  (cut) 

Velocity  Data  (m/sec) 

BBL  (1962) 

1 

McSkimin^ 

Obs. 

i 

1 

Calculated 

BBL  This  Work 

I  II 

T" 

1 

1 

Obs. 

Calc. 

0 

0 

J  a(Y) 

T 

1 

1 

6006 

6005 

0 

0 

b(Y) 

1 

1 

1 

4323 

4322 

0 

0 

c(Y) 

3892 

| 

3914 

3915 

3901 

1 

1 

3918 

3917 

0 

3  5°  1 5 1 

■  c(AT) 

3320 

1 

3320 

3322 

3246 

1 

1 

0 

45 

1  3 

1 

| 

7121 

7021 

0 

45 

1  b 

1 

1 

3981 

3981 

0 

45 

1  C 

1 

1 

3427 

3427 

0 

90 

|  a(Z) 

1 

1 

6319 

6319 

0 

90 

1  b,  c(Z) 

1 

1 

4687 

4688 

20 

34°20' 

1  oUT) 

3386 

1 

3566 

3551 

3467 

1 

30 

0 

1  a(X) 

5690 

1 

3746 

5746 

5643 

1 

5749 

5749 

30 

0 

1  b(X) 

1 

1 

5114 

5114 

30 

0 

1  c(X) 

1 

1 

I 

3298 

3298 

30 

10 

a 

5780 

1 

5960 

5956 

5830 

1 

30 

30 

b 

4096 

i 

4100 

4097 

4104 

I 

45 

34.  5° 

c(RT) 

4080 

i 

4094 

4080 

4117 

1 

45 

35 

1  b 

4416 

1 

4540 

4533 

4417 

1 

47.  5° 

33 

1  C 

4202 

1 

4176 

4162 

4205 

1 

50 

32 

1  C 

4212 

1 

4240 

4223 

4271 

1 

50 

33 

1  0 

4316 

1 

4200 

4183 

4228 

1 

50 

38 

1  b 

4694 

1 

4776 

4767 

4684 

1 

55 

47 

1  b 

5038 

1 

5030 

5021 

5004 

1 

60 

49°13 ' 

1  b(BT) 

5080 

1 

5060 

5072 

5089 

1 

Number  of  test  points 

14 

1 

11 

Least  Square  Variance 

7570 

7186 

2198 

1 

_1_ 

0.  5 

Table  8.  Elastic  Constants  of  Quartz 


pq 

1 

1 

c  in  units  of  109  N/m2 

pq  ' 

1  1 

McSkimin  ' 

BBL  (1962)  Data  1 

BBL  This  Work  1 

_ L 

BBL  (1961) 

n 

f  -  1 

|  8G.  80  | 

86.  74 

83.  63  1 

86.  3 

13 

!  11.91  1 

1  . 

11.  91 

-0.88  1 
| 

6.3 

14 

1  1 
|  -18.04  , 

-17.  91 

-18.88  | 

-17.  1 

33 

•  105.75  1 

1 

107.  2 

77.  60  1 

90.  3 

44 

1 

,  58.  20  1 

57.  94 

57.32 

66.  0 

6(5 

1  39.88  | 

39.  88 

39.  58  1 

32.  2 

ell 

=  0.  171  eH  -  -0.0406  (C/m2) 

E11 

=  39.  21  X  10'12  e 

33  *  41-03  X 

10‘12  (F/m) 

2 

ell 

/eu  ^  0.  77  X  109  (N/m2)  (McSkimin) 

P  = 

2650  kg/m3  (BBL) 

2648.  5  kg/m 

3  (McSkimin) 

McSkimin  et  al  have  applied  a  least  squares  minimization  to  the  experimental 

data,  and  the  derived  c  values  in  Table  8  are  the  result  of  this  analvsis.  They 
pq 

minimized  the  variance 

n 

a2  (1/n)  ^  G.  (v  -  v  .  )2  (41) 

/  -«  t  exp  calc  l 

il 

where  the  summation  is  over  the  number  of  experimental  points  n.  The  weighting 

2 

factor  G.  is  either  G.  =  1  or  G.  =  I/veXp-  If  the  experimental  points  are  equally 

reliable  on  an  absolute  basis,  then  G.  =  1  is  the  appropriate  choice.  If  the  exper- 

1  2 
imental  points  are  equally  valid  on  a  percentage  basis,  one  chooses  G.  =  I/veXp* 

For  this  data  set,  the  differences  between  the  two  approaches  are  insignificant. 


8.  McSkimin,  H.J.,  Andrcalch,  P.  ,  and  Thurstoir  R.  N.  (19(55)  Elastic  moduli 
of  quartz  versus  hydrostatic  pressure  at  25c'C'  and  -195.  8°C,  J.  Appl. 
Phys.  3(5:1524-1(532. 
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The  variances  for  the  various  data  sets  are  listed  at  the  bottom  of  Table  7. 

The  listed  velocity  values  are  rounded  off  to  the  nearest  integer,  whereas  the 

computations  were  carried  out  to  the  appropriate  number  of  decimal  places. 

Hence  there  are  small  differences  between  the  listed  variances  and  ones  computed 

from  the  listed  velocities.  The  close  agreement  between  calculated  and  experi- 

2 

mental  velocities  for  the  McSkimin  data,  as  reflected  in  a  =0.  5,  is  not  surpris- 

9 

ing.  McSkimin  used  10  of  the  11  test  points  to  derive  6  cpq  values  and  Cjj/ejj, 
and  one  does  expect  a  close  fit.  It  is  more  questionable  whether  7  fitting  param¬ 
eters  and  10  data  points  constitute  a  statistically  significant  set,  especially  if 
these  points  are  chosen  along  the  uncoupled  directions. 

The  variance  of  the  BBL  data  set,  using  BBL  (1962)  derived  material  con- 
2 

stants,  is  a  -  7186.  BBL  (1962)  have  not  performed  a  least  squares  analysis  to 
their  14  point  data  set,  and  from  their  publication  the  methodology  used  to  derive 
c^  is  unclear.  We  have  subjected  the  BBL  data  set  to  a  least  squares  minimiza¬ 
tion  routine,  and  the  results  of  the  best  fit  to  the  experimental  data  are  listed  in 

Table  7  under  "This  Work  II"  and  in  Table  8  under  "This  Work".  We  note  that  the 

2  2 

least  squares  fit  analysis  reduces  the  initial  a  -  7186  to  a  =  2198,  with  corre¬ 
sponding  differences  in  c  The  major  c^  changes  occur  in  Cjg  and  c^g,  the  two 
constants  not  accessible  to  uncoupled  mode  piezoelectric  excitation. 

We  have  performed  several  computational  efforts  trying  to  improve  the 
BBL  (1962)  data  fit.  One  included  the  assumption  that  one  of  the  14  test  points  is 
inaccurate.  We  systematically  neglected  one  test  point  at  a  time,  and  re-fitted 
the  data.  There  were  some  improvements,  but  none  of  them  dramatic.  The  best 
initial  13  point  data  fit  is  obtained  by  omitting  (30,  10)°,  with  the  variance  reduced 

to  ct2  -  5363,  an'4  the  best  final  fit  is  obtained  by  neglecting  the  AT-cut  (0.35.25)°. 
2 

w  ith  rj  -  1360.  However,  the  AT-cut  is  the  point  that  BBL  (1962)  assumes  to  be 

in  perfect  coincidence.  Another  supposition  was  that  BBL  (1962)  has  mislabeled 

one  of  the  modes.  The  mode  assignment  of  each  test  point  was  systematically 

2 

changed,  one  at  a  time,  and  the  data  re-fitted.  In  all  cases  the  initial  a  =  7186 
value  increased,  confirming  the  accuracy  of  the  BBL  (1962)  mode  assignments. 

There  is  another  BBL  data  set.  Ref.  6,  designated  as  BBL  (1961),  consisting 
of  17  rotations  with  47  test  points.  Some  test  points  of  the  two  sets  are  iden¬ 
tical,  and  some  are  listed  in  BBL  (1962)  with  slightly  different  9-angles.  Also, 
all  (4>.  9)  combinations  of  this  data  set  are  given  as  integers,  without  any  error 

bars,  a  highlv  unrealistic  situation.  Nevertheless,  we  have  performed  a  least 

2 

square  minimization  of  the  47  test  points  and  the  initial  o  =  5840  can  be  reduced 
2 

to  a  =  2903,  the  same  order  of  magnitude  as  the  BBL  (1962)  data  fit.  The  corre¬ 
sponding  c  values  for  this  fit  are  listed  in  Table  8  under  "BBL  (1961)". 
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7.  ANALYSIS  AND  CONCLUSIONS 


The  knowledge  of  the  elastic  constants  of  quartz  and  their  temperature 

coefficients  is  imperative  in  describing  the  frequency  and  temperature  response 

of  the  resonator  device.  In  addition  to  the  elastic  constants  c  ,  BBL  (1962) 

P9 

also  derived  first,  second,  and  third  order  temperature  coefficients  of  the  elastic 
constants,  Tn(Cpq),  n  =  1,  2,  and  3,  from  a  set  of  23-24  crystallographically 
doubly  rotated  resonators.  The  coefficients  were  obtained  by  fitting  the  static 
frequency-temperature  curves  in  the  -200  to  +200°C  range  by  the  cubic 
T^lfHT  -  T  )n  relationship,  where  Tn(f)  are  the  temperature  coefficients  of  fre¬ 
quency.  We  have  shown  that  there  are  serious  discrepancies  in  the  original 
analysis  of  the  BBL  (1962)  data  set,  and  this  raises  issues  regarding  the  applica¬ 
bility  of  currently  accepted  c  and  ^(c^)  values  to  crystallographically  doubly 
rotated  cuts. 

The  large  differences  in  c^^  are  especially  disturbing.  This  constant,  as 
indicated  in  Table  6,  is  the  easiest  to  determine  by  non-piezoelectric  means  from 
(</>,  90)°,  and  it  should  be  the  most  accurate.  Yet,  it  shows  the  largest  variation 
when  computed  from  an  arbitrary,  crystallographically  doubly  rotated,  piezo- 
electrically  excited  resonator  set.  At  the  same  time,  c^.,  computed  by  the  same 
technique  from  (<ti,  90)°,  does  agree  very  well. 

Our  major  concern  regarding  Tn^cpq)  relates  to  the  fact  that  the  computation 
of  the  first,  second,  and  third  order  temperature  coefficient  involves  a  knowledge 
of  the  elastic  constants,  and  that  the  higher  order  coefficients  also  depend  on  the 
values  of  the  lower  order  coefficients.  Errors  are  cumulative,  and  c  inaccu- 

pq 

racies  are  reflected  in  all  temperature  coefficients.  Similarly,  Tj(Cp^)  errors 

invalidate  T„(c  )  results,  and  these  in  turn  affect  T_(c  ).  The  BBL  (1962) 

2  pq  3  pq 

derived  temperature  coefficients  are  often  in  sharp  disagreement  with  previous 
results,  but  are  claimed  to  be  more  accurate  than  ones  determined  by  earlier 
investigators.  The  problems  associated  with  the  BBL  (1962)  data  set  regarding 
agreement  with  accepted  c  values  raises  more  serious  questions  with  respect 
to  the  BBL  (1962)  derived  temperature  coefficients. 

There  are  several  possibilities  which  may  account  for  the  BBL  (1962)  dis¬ 
crepancies.  Strictly  speaking,  the  derived  material  constants  are  not  necessarily 
the  intrinsic  values  and  their  applicability  to  other  situations  needs  to  be  estab¬ 
lished.  Geometric  factors  may  influence  results,  and  sets  of  values  derived  from 

one  configuration  and  experimental  procedure  may  not  automatically  be  applied  to 

g 

other  modes  of  operation.  Adams  et  al  also  derived  a  set  of  temperature 


9.  Adams,  C.A.,  Enslow,  G.  M. ,  Kusters,  J.A.,  and  Ward,  R.W.  (1970) 
Selected  topics  in  quartz  crystal  research,  24th  Annual  Symposium  on 
Frequency  Control,  AD  746210, 


30 


coefficients  of  elastic  constants,  and  they  caution  that  their  T  (c  )  data  is 

n  pq 

"exactly  applicable  to  crystals  having  the  same  physical  characteristics  as  the 
ones  in  this  study;  that  is,  same  diameter,  thickness,  contour,  etc."  BBL  (1963), 
Ref.  10  also  lists  factors  that  influence  the  temperature -frequency  behavior  of 
the  resonator. 

McSkimin's  experiments  were  carried  out  on  0.5-cubic-inch  quartz  samples, 
whereas  BBL  obtained  the  material  from  several  sources.  The  BBL  crystals  had 
large  thickness  to  diameter  ratio  variations,  and  included  plano-plano,  beveled 
at  the  edges,  and  plano-convex  geometries.  Sample  set  divergencies  extended  to 
fabrication  processes,  support  structure,  electrode  size,  air  gap  holder  and 
resonator  enclosure,  drive  level,  and  operations  at  different  overtones,  all  of 
which  affect  the  measured  frequencies.  Recent  discussions  with  A.  Ballato  tend 
to  support  the  argument  that  for  the  BBL  sample  set  they  lacked  an  accurate 
knowledge  of  (^,  9).  The  oriented  plates  and  finished  resonators  were  acquired 
from  outside  sources,  and  there  was  no  independent  check  on  (d.  9). 

Another  possibility  for  the  discrepancies  relates  to  the  fact  that  we  are  trying 
to  apply  the  theoretical  formulation  of  an  infinite  plate  to  a  finite  electroded  disk. 
This  is  very  simplistic  approach  that  works  well  when  symmetries  cancel  a  large 
number  of  terms,  but  does  not  describe  the  realities  of  doubly  rotated  disks.  It 
is  then  not  too  surprising  that  this  mixture  of  geometrical  configurations,  experi¬ 
mental  procedures,  and  simplistic  mathematical  model  yields  large  variances, 
and  deviations  from  predicted  results  can  always  be  attributed  to  some  of  these 
factors. 

We  are  confronted  with  two  issues:  (1)  the  validity  of  the  BBL  (1962)  deter¬ 
mined  c  owing  to  the  large  variance  associated  with  the  data  set,  and  (2)  that 

rM 

the  number  of  points,  14,  used  to  derive  the  least  square  fit  c  set  in  this  investi- 

PM 

gation,  is  still  statistically  insignificant.  Consequently,  despite  the  reduced 

variance,  there  is  no  assurance  that  this  c^  set  gives  a  more  realistic  prediction 

of  the  frequencies  of  doubly  rotated  cuts.  The  BBL  (1962)  c  values,  confirmed 

Pm 

by  independent  measurements  by  McSkimin,  are  in  general  agreement  with 
previous  data  accumulated  by  many  investigators  over  the  last  50-70  years. 
However,  all  previous  investigations  utilized  uncoupled  mode  orientations,  and. 
to  the  best  of  our  knowledge,  BBL  (1962)  was  the  first  group  trying  to  derive  c 

PM 

from  widely  rotated  ($,  9)  cuts.  Thus,  despite  all  the  reservations,  one  is  forced, 
at  the  present  time,  to  retain  and  accept  the  BBL  (1962),  or  the  McSkimin.  c 

PM 

values.  In  this  context,  the  least  squares  fit  derived  c  ,  listed  in  Table  8.  can 

pq 


10.  Bechmann,  R.,  Ballato,  A.D.,  and  Lukaszek,  T.J.  (1963)  Higher-Order 
Temperature  Coefficients  of  the  Elastic  Stiffnesses  and  Compliances  of 
o-Quartz,  USAELRDL  TR  2261. 


-  -3 


be  regarded  as  an  indication  that,  despite  all  previous  work  on  this  subject,  the 

derived  c  are  inconsistent, 
ptl 

This  investigation  points  out  the  need  for  the  experimental  redetermination 
of  the  elastic  constants  of  quartz  and  their  temperature  coefficients  on  geometric¬ 
ally  consistent  sample  sets,  on  a  statistically  significant  number  of  doubly  rotated 
crystals,  bulk  and  disk  configurations,  air  gap  and  resonator  enclosures,  and 
piezoelectric  and  other  means  of  excitations.  It  also  calls  for  detailed  documen¬ 
tation  of  sample  and  experimental  procedures.  A  new  experimental  investigation 
will  also  afford  the  opportunity  to  evaluate  novel  factors,  discovered  during  the 
last  twenty  years,  which  influence  resonator  behavior. 
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Appendix  A 

Transverse  Vibrations  Along  the  Xf-Axis 


7 

Slobodnik,  Delmonico,  and  Conway  define  the  acoustic  wave  propagation  in 
terms  of  a  fixed  coordinate  system  (x,  y,  z)  rotated  axes  (1,  2,  3),  Euler  rotation 
angles  (X,  p.  6),  and  a  phase  velocity  vector  that  always  propagates  along  the  1-axis. 
We  wish  to  specialize  the  Euler  angles  to  describe  a  doubly  rotated  Y-plate  with 
its  thickness  direction  normal  to  the  y-axis,  and  with  rotation  angles  (4>,6),  as 
depicted  in  Figure  1. 

The  first  Euler  angle  X  denotes  a  rotation  around  the  z-axis,  the  second  Euler 
angle  p  denotes  a  rotation  around  the  1-axis,  that  is,  the  rotated  x-axis,  and  the 
third  Euler  angle  (j  denotes  a  rotation  around  the  3 -axis,  that  is,  the  rotated  z- 
axis.  In  this  formulation,  the  initial  position  of  the  plate,  depicted  in  Figure  Al.a, 
is  that  of  an  X-plate,  thickness  direction  normal  to  the  crystallographic  x-axis. 

The  angle  X  =  90°  rotates  the  X-plate  into  a  Y-plate.  For  our  purposes, 

X  =  90  +  This  rotates  the  X-plate  past  the  Y-plate  position,  and  d  becomes  the 
rotation  angle  of  the  Y-plate  around  the  z-axls.  The  position  of  the  coordinate 
system  after  the  first  Euler  angle  rotation  is  depicted  in  Figure  A  l.b.  The 
second  Euler  .ingle  p  is  a  rotation  around  the  1-axis.  This  angle  has  to  be  chosen 
in  a  manner  that  will  align  the  3-axis  with  the  x-y  plane.  This  is  necessary  since 
the  third  Euler  angle  0,  rotation  around  the  3 -axis,  must  be  identical  with  the 
x'-axis  rotation  depicted  in  Figure  1,  and  the  x'-axis  lies  in  the  x-y  plane,  p  is 
then  chosen  as  p  =  90°,  and  the  position  of  the  coordinate  system  after  the  second 
Euler  angle  rotation  is  depicted  in  Figure  Al.c.  The  position  of  the  coordinate 
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system  after  the  third  Euler  angle  rotation,  a  rotation  0  around  the  3 -axis,  is 
depicted  in  Figure  Al.d. 

The  Euler  angles  for  a  doubly  rotated  Y -plate  are  then  given  by 
(X,p,0)  (90  f  90,  6) 

For  these  angles  the  Euler  transformation  matrix  reduces  to 


-sin  <t>  cos  0  cos  <t>  cos  0  sin  0 
V...  =  sin  <t>  sin  0  -cos  <t>  sin  0  cos  0 


The  rotated  material  properties  e!.,  e!  ,  and  c'  are  then  determined  from 

t.T  ip  pq 

Eqs.  (19)-(21)  and  transformation  matrix  of  Eq.  (A2).  Slobodnik  has  not  derived 
explicit  algebraic  expressions  for  the  rotated  properties,  but  numerical  results 
can  be  easily  generated  by  computerized  DO  LOOP  techniques.  It  can  be  shown 
that  there  is  equivalence  between  the  expressions  listed  in  Table  3  and  the  ones 
derived  from  Eq.  (A2).  Table  A1  lists  the  notations  of  this  equivalence. 
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Table  Al.  Comparison  of  " Rotated" 
Material  Constants 


This  Work 

Slobodnik 

c22 

CU 

c24 

C  * 
c16 

c26 

c  * 
c  15 

c44 

c  * 
c  66 

c  * 

C4G 

f*  • 

56 

c66 

C55 

e22 

e'll 

e24 

e16 

e26 

e'l5 

E  22 

£ii 

For  propagation  along  the  x^axis,  the  eigenvalue  determinant,  Eq.  (11),  is 
expanded  with  v  ~  1 ,  and  one  obtains 


ril'r  ri«  ri5 

rKi  '  r  r5(i 

ri5  r5fi  r55‘r 

-  0 

(A3) 

The  coefficients  of  the  eigenvalue  equation 

r3  +  A2r2  +  Ajr  +  a0  =  o 

(A4) 

become 

A2  =  '<ri  1  +  r55  +  r66) 

(A  5) 

A1  (riir55  ’  r  15'  +  (riir6G 

-  T  l6)  +  <r55r66  ■  r56* 

(A  6) 

A  -  r.  ,r?r  +  r,cr?r  +  r„rr?c  -  r.  .r,,.rcc  -  2r1crlcr,c 

0  11  56  55  1G  6 B  15  11  55  66  15  16  56 

(A  7) 
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A1  show  that,  in  terms  of  rotation 
o  the  crystallographic  axes,  both  formu- 
lts  for  the  cubic  equation,  and  the  eigen- 


